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Abstract. Series expansion methods are used to study directed bond percolation 
clusters on the square lattice whose lateral growth is restricted by a wall parallel to 
the growth direction. The percolation threshold Pc is found to be the same as that for 
the bulk. However the values of the critical exponents for the percolation probability 
and mean cluster size are quite different from those for the bulk and are estimated by 
j3i = 0.7338 ± 0.0001 and 71 = 1.8207 ± 0.0004 respectively. On the other hand the 
exponent Ai = /3i + 71 characterising the scale of the cluster size distribution is found 
to be unchanged by the presence of the wall. 

The parallel connectedness length, which is the scale for the cluster length 
distribution, has an exponent which we estimate to be = 1.7337 ± 0.0004 and 
is also unchanged. The exponent ti of the mean cluster length is related to Pi and 
z^iii by the scaling relation i>m — Pi + ti and using the above estimates yields ti — 1 
to within the accuracy of our results. We conjecture that this value of ti is exact and 
further support for the conjecture is provided by the direct series expansion estimate 
Ti = 1.0002 ± 0.0003. 
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Recently exact results have been obtained for directed compact clusters on the 
square lattice near a wall [0, ^ ^. Such clusters are similar to ordinary percolation 
clusters except that they cannot branch and have no holes. These simplifying features 
allow several of the usual percolation functions to be derived analytically and the 
corresponding critical exponents have integer values. One of the main conclusions from 
these results was that although the moments of the cluster size and length distributions 
have exponents which change on introducing a wall parallel to the growth direction the 
exponents for the size and length scales remain the same. The other was that growth 
parallel to the wall is rather special in that any bias away from the wall results in bulk 
exponents. Similarly any bias towards the wall leads to wet wall exponents 

In this paper we find that the first of these conclusions extends to directed bond 
percolation. The exponents for directed percolation are not known exactly but numerical 
results show that, even in the absence of a wall, they are generally far from being integer 
and there is some doubt as to whether they even have rational values 0. An interesting 
possibility raised by our results is that the mean cluster length in the presence of a 
wall parallel to the growth direction is an exceptional case and has the integer exponent 
Ti = 1. Direct evidence for this value is provided by our analysis of the low density 
series expansion for the mean cluster length. Further support is provided by the scaling 
relation 

Pi + ri = z^iii (1) 

together with series expansion estimates of (3i and i/iy. Here the subscript 1 on an 
exponent indicates its value in the presence of a wall. This relation is less well known 
than the one for the cluster size distribution, namely 

Pi + 71 = Ai (2) 

and is derived below. First we define the model and introduce some notation. 

The directed square lattice may be described as having sites which are the points in 
the t — x plane with integer co-ordinates such that t > and t + x is even. There are two 
bonds leading from the general site (t, x) which terminate at the sites {t + l,x±l). All 
bonds have probability p of being open to the passage of fluid and the source is placed 
at (0,0). This will be known as the bulk problem. A wall will be said to be present 
if the bonds leading to sites with x < are always closed. The probability that fluid 
reaches column t but no further will be denoted by rt{p) and in this event the origin 
will be said to belong to a cluster of length t. 

The percolation probability, the probability that the origin belongs to a cluster of 
infinite length, is defined by 

oo oo 

P{P) = 1-11 ^t(p) = - nip)) ~ (p - Pcf for p p+. (3) 

t=0 t=0 
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If we suppose that the length distribution has the scaling form 

(4) 

then if 

Up)--\Pc-pn (5) 

substitution in @ yields 

a = 1 + — . (6) 

Z/|| 

The mean cluster length is defined by 

oo 

Tip) = Y.tnip) (7) 

t=0 

and using (H) we find that 

T(p)~bc-p|-" (8) 



where 



z/||-/?. (9) 



The same argument holds in the presence of the surface and leads to (|1|). There is a close 
correspondence between the above derivation and that of (0) given in 0. To obtain (|^) 
it is only necessary to replace rt{p) by the cluster size distribution Ps{p), ^\\{p) by the 
scaling size cr(p), which diverges with critical exponent A, and T{p) by the mean cluster 
size S{p) which diverges with exponent 7. 

The mean size and the parallel and perpendicular scaling lengths are obtained from 
the pair connectedness function C{t,x;p) which is the probablity that there is an open 
path from the origin to the site {t,x). The moments are defined by 

/^m,n(p)= E^"^"C'(t,X;p) (10) 

sites 

in terms of which S{p) = /ioo(p)- Assuming a scaling form for C{t,x;p) similar to 
where x is scaled by ^±{p), it follows that 

^Wip) TT and C±{P) f^- (11) 

The series expansions are obtained by a transfer matrix method similar to that used for 
the bulk lattice 0] and the details of the implementation in the presence of a wall will 
be given in a forthcoming paper |^. The state of column t is a specification of which 
sites in that column are wet and which are dry and the probability that state i occurs 
is denoted by Hi(t,p). The state in which all sites are dry is labelled i = 0. Essentially 
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the state vector of a given column is completely determined by that of the previous 
column and only one state vector need be held in the computer at any stage. C(t,x;p) 
is determined by summing 'Ki{t^p) over all states for which the site with co-ordinate x 
is wet and 

=7ro(t + l,p)-7ro(t,p). (12) 

Low density expansions in powers of p are obtained by noting that TT{t,p) = 0{p^) 
so that all of the above functions may be obtained to this order by computing the 
state vectors up to column t. We were able to derive the series directly up to a maximal 
column tm = 49. However, these series can be extended significantly via an extrapolation 
method similar to that of As an example, consider the series for the average cluster 
length T(p). For each t < we calculate the polynomials Tj(p) = J2t'=o'^'^t'iP) correct 
to 0{p™). As already noted these polynomials agree with the series for T{p) to 0{p^). 
Next, we look at the sequences df^s obtained from the difference between successive 
polynomials 

Tt+i{p)-Z{p)=p'+'Y.dt,sP'. (13) 

s>0 

The first of these correction terms dtfl is often a simple sequence which one can readily 
identify. In this case we find the sequence 

-dtfi = 1, 2, 3, 6, 10, 20, 35, 70, 126, 252, 462, 924, . . . 
from which we conjecture 

d2tfl = 2d2t-i,o, f/2i-i = 475(t+ 1,-1/2), (14) 

where B{x,y) is the Beta function. The formula for dt^ holds for all the trn.~l values 
that we calculated and we are very confident that it is correct for all values of t. As 
was the case in p|] the higher-order correction terms dt^s can be expressed as rational 
functions of dtfl, 

''/'I f t-s\ 

dt,s = X! 7 {C'S,kdt~s+lfl + bs,kdt-s+2fl) + 5Z (^s,kdt-s+k+lfl- (15) 

k=l \ ^ / k=0 

From this equation we were able to find formulas for all correction terms up to s = 17 
and using T^g^p) we could extend the series for T{p) to 0{p^'^). A similar procedure 
allowed us to extend the series for S{p) and the parallel moments /ii,o(p) and fi2fl{p) 
to 0{p^'^), while the series for the first and second perpendicular moments, ^o,i{p) and 
f^o,2{p), were extended to 0{p^^). The resulting series are listed in table More details 
of the extrapolation procedure including the formulae for the various correction terms 
will appear in a later paper [^. 
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The only high density expansion we consider is that for the percolation probability 
which can be obtained from (|12|) and (|^) by noting that rt{p) = 0{q^) where q = 1 —p 
and k is the least integer > ^(t + 2). Thus for a given value of t the number of 
terms obtainable in the high density expansion is only about half as many as in the low 
density expansion. However, for computational purposes it is more efficient to derive the 
series expansion for P{q) directly via a transfer matrix technique. For the percolation 
probability we derived the series directly to Olq"^^) and obtained another 8 terms from 
the extrapolation procedure. The resulting series is listed in table ^ 

It is found from unbiassed approximants that the estimates of pc agree with the 
bulk value @], Pc = 0.6447002 ± 0.0000005 obtained from longer series and we therefore 
bias our exponent estimates with this value. This value of pc was obtained from low 
density series and is a refinement of that obtained from analysis of the shorter series for 
P{q) ^ which gave pc = 0.6447006 ±0.0000010. Data obtained from T(p), the parallel 
moments and P{q) is shown in tables ^, Hand|^. The exponent of /ioo(p) was estimated 
from the series for {S{p) — l)/p which is the mean size of the cluster connected to the 
site (1,1); this gave better convergence. We have also analysed the first and second 
perpendicular moment of the pair connectedness and series for C,\\{p) and C,±{p) obtained 
from (|ll]) using the first and second moments. In the analysis of P{q) we used standard 
DLog Fade approximants while the remaining series were analysed using first and second 



order inhomogeneous differential approximant |11 



In table ^ the columns headed L = result from the standard DLog Fade analysis 
and give ri = 1 to three decimal places although most of the entries are slightly above. 
This conclusion is not altered by looking at inhomogeneous approximants (the first few 
of which we have included in table |^) or second order approximants. Using the slightly 
smaller value pc = 0.6446980 gave the better converged result ti = 1.00004 ± 0.00004. 

We turn now to the indirect evidence for ti = 1 via the scaling relation (|l]). The 
value z/i|| = 1.7337 ± 0.0004 was obtained by analysing the series for /i2,o(p)//Ui,o(p) and 
is consistent with the value obtained by subtracting the value of the exponent of /ii,o 
from that of /i2,o- It is clearly equal to the corresponding bulk exponent, as in the case of 
compact percolation, and we use the more accurate bulk estimate in deriving ti below. 
The corrections to scaling in the case of the percolation probability appear to be very 
close to analytic, and the standard Fade estimate of /3i (table H) should be accurate. 
Combining the values of z/y and jSi gives ri = 1.0000 ± 0.0002 which agrees with the 
direct estimate. 

Other exponent values obtained from the analysis of various series are collected 
together in table ^ where previous estimates for the bulk problem and exact results 
for compact percolation are also given. As usual the error bars are a measure of the 
consistency of the higher order approximants and are not strict bounds. The estimate 
(3 = 0.27643 ±0.00010 of |TD| has been adjusted slightly upwards to allow for the change 
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in pc- In estimating the exponents we rely both on the analysis of the series yielding 
a particular exponent and estimates obtained using scaling relations. In some cases we 
also use the more accurate bulk exponent estimates. A case in point is the exponent 
7i. From the Dlog Fade approximants in table ^ one would say that the direct estimate 
from the series for {S{p) — l)/p favours a value of 71 ~ 1.8211 with a rather large spread 
among the approximants. However, the better converged estimates of 7i + 2i/i|| ~ 5.2881 
together with the bulk estimate of i>\\ leads to 71 ~ 1.8205. In this case second order 
differential approximants to S{p) are better converged and favour 71 ~ 1.8207. Taking 
all the evidence into account including our belief that Ai takes on the bulk value we 
arrived at the estimate for 71 quoted in table The estimate of r is derived from the 
scaling relation (H). Analysis of the bulk expansions @, ||, |10] showed that corrections 
to scaling were close to analytic, as they are here. 

The values of i'i± and Ai (obtained from the scaling relation @), as well as z/iy, 
are clearly the same as those for the bulk. The scaling size and both scaling lengths are 
therefore unchanged by the introduction of the wall. We also note that the hyperscaling 
relation, with D the dimension of space perpendicular to the preferred direction t {=1 
for the square lattice), 

+Du^ = (3 + A, (16) 

which is satisfied by the bulk exponents apparently fails on the introduction of a wall. 



We now consider the possibility of rational exponents. As previously noted |T0 
there is no simple rational fraction whose decimal expansion agrees with the estimate 
of p. The same is true for other exponent estimates in table |^. In particular we note 
that our estimates of the bulk exponents z/|| and z/^ differ by 0.03% from the rational 
fractions u\\ = 26/15 = 1.733 333..., and u± = 79/72 = 1.097 222... suggested by 
Essam et al 0. We believe this to be a significant difference given the high precision 
of our results. However, the suggested rational fraction 7 = 41/18 = 2.277 777. . . and 
the value of A = 613/240 = 2.554 1666 . . ., which follows from the above rational values 
by scaling, are generally still within our estimated the error bounds. The fraction for 
A is not very appealing though and assuming that both exponents have these values 
then scaling implies the even less convincing result (3 = 199/720 = 0.276 388 . . . which 
is however just consistent with our estimated value. 

If we assume that ri = 1 is exact and that the values of A, z/y and i>± are the same 
with and without a wall then all of the other surface exponents are determined by scaling 
together with the values of any three bulk exponents. The surface exponents calculated 
in this way are presented in table |^ for comparison with the estimated values of table 
^ as a measure of the overall consistency of our results. The bulk exponents used were 
7 = 41/18 and the bulk estimates of and z/_l. Excellent agreement is observed. 

Our findings may be summarised as follows. Firstly we have found that the scaling 



7 



size and both scaling length exponents are unchanged by the introduction of a wall 
parallel to the preferred direction. Also we have examined the widely held view that two 
dimensional systems should have rational exponents. The high precision data presented 
here is consistent with the results Ti = 1 and 7 = 41/18. However there are no such 
simple fractions which are in agreement with our estimates of z/|| and z/_l. Given that 
directed percolation is not conformally invariant, and that the expectation of exponent 
rationality is a consequence of conformal invariance, this is perhaps not surprising. 
The precise numerical work reported and quoted in this paper therefore supports the 
conclusion that the critical exponents for non translationally invariant models should 
not in general be expected to be simple rational numbers. The cluster length exponent 
Ti and the exponent 7 appear to be exceptional cases. 
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Tables and table captions 



Table 1. Low density expansions in powers of p, row n is the coefficient of p". 
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Table 2. High density expansion for the percolation probabihty P{q) = J2^nQ"- 
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Table 3. Differential approximant analysis of the mean length scries. The table 
shows biased first order inhomogeneous approximant estimates of ri . L is the degree 
of the inhomogeneous polynomial. For i = the entries are from biased Dlog Fade 
approximants. 
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1 


.00028 
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1 


.00023 


32 


1 


.00023 


1 00022 


1 00022 


1 


.00023 


1 


00022 


1 


.00022 


33 


1 


.00022 


1 00022 




1 


.00022 
















L = 2 










L = 3 






N 


[N-1,N] 


[N,N] 


[N+1,N] 


[N-1,N] 




[N,N] 


[N+1,N] 


22 


1 


.00010 


1.00001 


0.99954 


1 


.00010 





.99996 


1 


.00014 


23 





.99962 


1.00004 


1.00009 


1 


.00014 


1 


.00008 


1 


.00011 


24 


1 


.00009 


0.99996 


1.00014 


1 


.00012 


1 


.00017 


1 


.00015 


25 


1 


.00014 


1.00017 


1.00144 


1 


.00015 


1 


.00019 


1 


.00034 


26 





.99986 


1.00028 


1.00030 


1 


.00038 


1 


.00030 


1 


.00029 


27 


1 


.00030 


1.00028 


0.99995 


1 


.00030 


1 


.00029 


1 


.00018 


28 


1 


.00014 


1.00020 


1.00020 


1 


.00019 


1 


.00020 


1 


.00022 


29 


1 


.00020 


1.00020 


1.00025 


1 


.00023 


1 


.00022 


1 


.00022 


30 


1 


.00030 


1.00022 


1.00023 


1 


.00022 


1 


.00023 


1 


.00023 


31 


1 


.00023 


1.00023 


1.00022 


1 


.00023 


1 


.00022 


1 


.00023 


32 


1, 


.00023 


1.00023 




1 


.00023 
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Table 4. DLog Pade analysis of the moments of the pair connectedness. The table 
shows biased approximant estimates of the critical exponents of the moments i^oo{p), 
Miob) and /U2o(p). 

71 7i+z^i|| 7i + 2i^i|| 



N [N-1,N] [N,N] [N+1,N] [N-1,N] [N,N] [N+1,N] [N-1,N] [N,N] [N+1,N] 



22 


1, 


.82381 


1 


.82760 


1 


.81953 


3 


.55492 


3 


.55555 


3, 


.55458 


5, 


.28807 


5, 


.28807 


5, 


.28808 


23 


1, 


.82010 


1 


.82593 


1 


.82355 


3 


.55466 


3 


.55478 


3, 


.55472 


5, 


.28809 


5, 


.28769 


5, 


.28807 


24 


1, 


.82364 


1 


.82094 


1 


.71766 


3 


.55473 


3, 


.55479 


3, 


.55462 


5, 


.28808 


5, 


.28804 


5, 


.28809 


25 


1, 


.77437 


1 


.81558 


1 


.82399 


3, 


.55466 


3. 


.55456 


3, 


.55457 


5, 


.28809 


5, 


.28809 


5, 


.28809 


26 


1, 


.82511 


1 


.82793 


1 


.82424 


3 


.55457 


3 


.55456 


3, 


.55459 


5, 


.28805 


5, 


.28809 


5, 


.28809 


27 


1, 


.82524 


1 


.82063 


1 


.82122 


3 


.55460 


3 


.55459 


3, 


.55459 


5, 


.28809 


5, 


.28808 


5, 


.28806 


28 


1, 


.82124 


1 


.82097 


1 


.82078 


3, 


.55460 


3, 


.55459 


3. 


.55459 


5, 


.28807 


5, 


.28809 


5. 


.28516 


29 


1, 


.82079 


1 


.82090 


1 


.82396 


3 


.55460 


3 


.55453 


3, 


.55458 


5, 


.28804 


5, 


.28802 


5, 


.28819 


30 


1, 


.81878 


1 


.82098 


1 


.82107 


3, 


.55459 


3, 


.55458 


3, 


.55459 


5, 


.28804 


5, 


.28806 


5, 


.28799 


31 


1, 


.82108 


1 


.82108 


1 


.82107 


3 


.55473 


3 


.55454 


3, 


.55463 


5, 


.28805 


5, 


.28761 


5, 


.28799 


32 


1, 


.82108 


1 


.82104 


1 


.82106 


3 


.55425 


3 


.55471 


3, 


.55460 


5, 


.28805 


5, 


.28806 


5, 


.28808 


33 


1, 


.82106 


1 


.82104 






3, 


.55589 


3 


.55471 






5, 


.28803 


5, 


.28806 







Table 5. DLog Pade analysis of the percolation probability series. The table shows 
biased approximant estimates of /3i. 



N 


[N-1,N] 


[N,N] 


[N+1,N] 


8 


0.73406 


0.73406 


0.73408 


9 


0.73409 


0.73409 


0.73408 


10 


0.73409 


0.73403 


0.73369 


11 


0.73389 


0.73388 


0.73385 


12 


0.73389 


0.73381 


0.73382 


13 


0.73382 


0.73381 


0.73382 


14 


0.73383 


0.73381 


0.73382 


15 


0.73382 


0.73382 


0.73379 


16 


0.73380 


0.73382 
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Table 6. Exponent values for compact and bond percolation. The bulk values for 
bond percolation are from except for (3 which is from adjusted for a small 
change in pc- The compact percolation results are from Q and references therein. 
Values in brackets are obtained from scaling formulae. The "with wall" value of 7 is 
from second order differential approximants. 





bond percolation 


compact percolation 


exponent 


with wall 


bulk 


with wall 


bulk 


T 


1.0002±0.0003 


(1.4573±0.0002) 





1 




0.7338±0.0001 


0.27647±0.00010 


2 


1 


7 


1.8207±0.0004 


2.2777±0.0001 


1 


2 




3.5546±0.0002 


4.0113±0.0003 


(3) 


(4) 


1 + 2i^|| 


5.2881±0.0002 


5.7453±0.0004 


(5) 


(6) 




1.7337±0.0004 


1.7338±0.0001 


(2) 


2 




4.014±0.002 


4.4714±0.0004 




(3) 


v± 


1.0968±0.0003 


1.0969±0.0001 




1 


A 


(2.5545±0.0005) 


(2.5542±0.0002) 


3 


3 



Table 7. Scaling values of the exponents for bond percolation calculated using ri = 1, 
7 = and the bulk estimates of and v±. 



exponent with wall bulk 



r 




1 


1.4573 






0.7338 


0.27646 


7 




1.8204 


2.2778 


74 


-^11 


3.5542 


4.0116 


74 


- 2v\\ 


5.2880 


5.7454 






1.7338 


1.7338 


74 


- 2iy± 


4.0142 


4.4716 






1.0969 


1.0969 


A 




2.5542 


2.5542 



